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PREDICTION OF THE PROTON-TO-TOTAL TURBULENT HEATING IN THE SOLAR WIND 
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ABSTRACT 

This paper employs a recent turbulent heating prescription to predict the ratio of proton-to-total 
heating due to the kinetic dissipation of Alfvenic turbulence as a function of heliocentric distance. 
Comparing to a recent empirical estimate for this turbulent heating ratio in the high-speed solar 
wind, the prediction shows good agreement with the empirical estimate for R > 0.8 AU, but predicts 
less ion heating than the empirical estimate at smaller heliocentric radii. At these smaller radii, 
the turbulent heating prescription, calculated in the gyrokinetic limit, fails because the turbulent 
cascade is predicted to reach the proton cyclotron frequency before Landau damping terminates the 
cascade. These findings suggest that the turbulent cascade can reach the proton cyclotron frequency 
at R < 0.8 AU, leading to a higher level of proton heating than predicted by the turbulent heating 
prescription in the gyrokinetic limit. At larger heliocentric radii, R > 0.8 AU, this turbulent heating 
prescription contains all of the necessary physical mechanisms needed to reproduce the empirically 
estimated proton-to-total heating ratio. 
Subject headings: turbulence — solar wind 



1. INTRODUCTION 

In the development of a thermodynamic model of the 
heliosphere, a crucial issue is the identification of the var- 
ious physical mechanisms that play a role in determining 
the temperature profiles of the heliospheric plasma ions 
and electrons. Some of the effects influencing the mea- 
sured temperature profiles in the inner heliosphere are 
energy conservation in the spherically expanding plasma, 
heat conduction, Coulomb collisions, and plasma heating 
through the dissipation of solar wind turbulenceQ The 
incorporation of the first three of these effects is relatively 
well understood, but the heating of the heliospheric ions 
and electrons resulting from the dissipation of solar wind 
turbulence remains an area of active research. 

Although a number of studies have addressed the issue 
of plasma heat i ng by the dissipation of turbulence (see 
iCranmer et~ai1 (|2009D for a review of previous efforts), 
with the exception of a pioneering series of papers by 
Quataert and Gruzinov (lQuataertl l l998t IGruzinovlll998l : 
iQuataert fc Gruzinovill9991 ). the turbulent heating of the 
ions and electrons separately has only recently become a 
focus of interest. We review here three studies that have 
made progress in the investigation of turbulent plasma 
heating and the effect of this heating on the solar wind 
p roton and elec t ron te mperature profiles. 

iBreech et al.l (|2009f ) coupled an e xisting turbulence 
transport model for the solar wind dZank et al.l Il996t 
Matt haeus et all 119991: iSmith et al.l 120011: IBreech et al l 



20081 ) to separate radial evolution equations for the pro- 



ton and electron temperatures. The temperature equa- 
tions included the effects of spherical expansion, parallel 
electron heat flux, collisional energy exchange between 
protons and electrons, and separate turbulent heating 
rates for the ions and electrons. The kinetic dissipa- 
tion mechanisms that determine the partitioning of tur- 

1 Pickup ions are believed to significantly affect the heliospheric 
energy balance on ly in the outer heliosphere at R > 10 AU 
HBreech et al.ll200g| y 



bulent heating between protons and electrons were not 
addressed in this model, so the fraction of proton-to- 
total turbulent heating was set to a constant value, 
fp = Qp/{Qp + Qe)- The model achieved a reason- 
able accounting for the temperature data measured by 
Ulysses when the effect of Coulomb collisions was weak, 
the electron heat conduction was provided by an em- 
pirically determined function, and the partitioning of 
turbulent heating was taken as f p = 0.6. Note that 
the value of f p = 0.6 is consistent with other estimates 
of the fraction of p roton heating (jLeamon "etaTI IT999I : 
iStawarz et"aLll2009f ). 

In a complementary work, ICranmer et ail (|2009) con- 
structed a model following the internal energy evolution 
of the protons and electrons in the solar wind in an at- 
tempt to constrain empirically the required proton and 
electron turbulent heating in the solar wind. Assuming 
isotropic Maxwellian proton and electron velocity dis- 
tributions with shared bulk velocity, separate equations 
for the conservation of energy in a spherically expand- 
ing flow for protons and electrons were constructed to 
incorporate the effects of energy exchange by Coulomb 
collisions, parallel electron heat conduction, and turbu- 
lent heating of protons and electrons. Using analytic 
fits to Helios and Ulysses measurements for the proton 
temperature, electron temperature, and parallel electron 
heat conduction flux, the equations were solved for the 
turbulent proton and electron heating rates. The results 
were then combined to obtain an empirical estimate of 
the proton-to-total heating rate Q P /(Q P + Q e ) as a func- 
tion of heliocentric radius, reproduced as the dashed line 
in Figure [T] 

Based on a theoretical model of the turb ulent cascade 
of ene r gy in a weakl y collisional plasma (|Howes et al.1 
I2008al) . IHowesI (|2010f ) predicted the proton-to-electron 
heating ratio Q p /Q e resulting from the dissipation of 
Alfvenic turbulence. The key result of this study 
was an analytical prescription for the heating ratio 
Qp/Q e (/3p,T p /T e ), a function of only two plasma pa- 
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rameters, the proton plasma beta /3 P and the proton-to- 
electron temperature ratio T p /T e . The limits of validity 
of this heating ratio prediction were given as a constraint 
on the minimum scale of turbulent energy injection (as- 
suming an isotropic driving mechanism for the turbu- 
lence) . 

Thi s paper describes the application of the iHowesI 
( 2010) heating prescription to predict the proton-to-total 
turbulent heating rate Q P /(Q P + Q e ) for the high-speed 
solar wind and comp ares the resulting pred iction to the 
empirical estimate of ICranmer et all (|2009Q . 

2. PREDICTION OF TURBULENT HEATING 

This section describes the prediction of the ratio of the 
proton-to-total turbulent heating Q P /(Q P + Q e ) in the 
high-speed solar wind as a function of heliocen tric ra- 
dius R using the turbulent heating prescription bv lHowe"sl 
(2010). Since this prescription depends on the plasma 
parameters f3 p and T p /T e , it is necessary to construct a 
model of the high-speed solar wind to determine the vari- 
ation of these plasma parameters with heliocentric ra- 
dius. In ij'2.11 we describe the solar wind model. In £12.21 
we review the theoretical framework of low-frequency, 
anisotropic Alfvenic turbulence in a magnetized, weakly 
collisional plasma that underlies the turbulent heating 
prescription. This prescription is employed to predict 
the proton-to-total turbulent heating Q P /{Q P + Q e ) in 
In ij2.4[ we estimate the evolution of the width of 
the inertial range in the high-speed solar wind in order to 
verify the validity of the turbulent heating prescription 
in §tE\ 

2.1. Solar Wind Model 

We adopt the s ame specific model fo r the high-speed 
solar wind used by ICranmer et al.l (|2009f ) to facilitate the 
comparison to their empirical turbulent heating estimate. 
This model is used to specify, as a function of heliocentric 
radius R, the two key plasma parameters required by the 
turbulent heating prescription: the proton plasma beta 
f3 p and the proton-to-electron temperature ratio T p /T e . 

Analytic fits to in situ measurements of the high-speed 
solar wind (faster than 600 km s _1 ) from the Helios and 
Ulysses spac ecraft over the range 0.29 AU < R < 5.4 AU 
were used bv lCranmer et al.l (|2009f ) to generate equations 
for the proton and electron temperatures as function of 
heliocentric radius R, 

m ( -^r77 I = 0.9711 - 0.7988a; + 0.07062a; 2 (1) 
\ 10 5 K / w 

In (j^-^j = 0.03460 - 0.4333a: + 0.08383a; 2 , (2) 

where x = ln(i?/l AU). 

In addition to the proton temperature, we need to 
specify the form of the proton density and magnetic 
field strength to determine the proton plasma beta /3 p = 
&Tm p T p /B 2 . Following ICranmer et al.l (|2009fl . we take a 
proton density of the form 

n p {R) = n {R/l AU)" 2 (3) 

where no — 2.5 cm~ 3 . T he empirical turbulent heating 
constraints calculated by ICranmer et al.l (|2009l ) used a 



colatitude 9 — 15° to model the high-latitude Ulysses 
measurements. For the heliocentric distances covered by 
this model, the winding of the magnetic field into the 
Parker spiral for the high-speed streams at this colatitude 
is relatively weak, so a simple monopolar model for the 
magnetic field strength is a reasonable approximation, 

B{R) = B (R/1 AU)" 2 , (4) 

with B = 2.5 x 10~ 5 G. 

Using these functions for T p , T e , n p , and B in high- 
speed solar wind streams, we find that the proton plasma 
beta varies from /3 p = 0.92 at 0.29 AU to /3 p = 34 
at 5.4 AU, and the proton-to-electron temperature ra- 
tio varies from T p /T e = 3.9 at 0.29 AU to T p /T e = 1.3 
at 5.4 AU. 

2.2. Turbulent Cascade Model 

The heating prescription presented in IHowesI (|2010l ) 
is determined using a model for the turbulent cascade 
of energy in a ma gnetized, weakly collisional plasma 
(jHowes et al.ll2008al ). The cascade model determines the 
steady state form of the magnetic energy spectrum of 
Alfvenic fluctuations, based on three primary assump- 
tions: (1) the Kolmogorov hypothesis tha t the energy 
casca de is determined by local interactions (jKolmogorovi 
119411 ): (2) the turb ulence maintains a state o f critical 
balance at all scales (jGoldreich fc Sri dhar 1995|); and (3) 
the linear kinetic damping rates are applicable in the 
nonlinearly turbulent plasma. 

The dependence of the nonlinear energy transfer rate 
on the local turbule nt fluctuations in the cascade model 
(jHowes et al. 2008a) is inspired by the following theoret- 
ical picture of low- frequency, anisotropic Alfvenic turbu- 
lence in a magnetized, weakly collisional plasma ([Howesl 
[200l ISchekochihin et al J [2001 . The energy of Alfvenic 
fluctuations is injected into the turbulence isotropically 
at a scale much larger than the ion Larmor radius, 
Lo 3> pi , corresponding to an isotropic driving wavenum- 
ber kopi <C 1. Since the damping of Alfvenic fluctua- 
tions at this large scale by wave-particle interactions in a 
weakly collisional plasma is negligible, the turbulent fluc- 
tuations rise to sufficient amplitudes that nonlinear inter- 
actions between counter-propagating Alfven wave pack- 
ets transfer the turbulent fluctuation energy to smaller 
scales. This sets up a critically-balanced, anisotropic cas- 
cade of MHD Alfven waves over all scales down to the 
perpendicular scale of the i on Larmor rad ius, k±pi < 1 
(jGoldreich fc Sridhan [19951: IBoldvrevll27)05P l. Even in a 
weakly collisional plasma, the dynamics of this Alfven 
wave cascade i s rigorously described by the equations of 
reduced MHD (jSchekochihin et al.ll2009j) . At the perpen- 
dicular scale of the ion Larmor radius k±pi ~ 1, the tur- 
bulence transitions to a critically balanced, anisotropic 
cascade of kinetic Alfven waves over the perpendicular 
scales k±pi > 1. 

The range of scales traversed by the MHD Alfven 
wave cascade, between the driving scale and the ion 
Larmor radiu s scale, is common ly designated the "in- 
ertial range" (jKolmog orovi Il94"lh of MHD turbulence — 
i.e., the range of scales over which the effects of driv- 
ing and dissipation are negligible. For a sufficiently 

2 Perpendicular and parallel are defined with respect to the di- 
rection of the local mean magnetic field. 
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large inertial range, the anisotropy of the energy cas- 
cade (fc|| oc k 2 / 3 in the Goldreich-Sridhar theory, or 

1 /2 

fcii oc in the Boldyrev theory) leads to turbulent 

fluctuations, at the transition to the kinetic Alfven wave 
cascade at k±pi ~ f, that are highly elongated along 
the direction of the local magnetic held, k»/k± <C 1. 
Such anisotropic fluctuations are optimally described by 
a low- fr equency expansion of kinetic th e ory called gyroki- 
netic s (iRutherford fe Friemanl [l968t iFrieman fe Chen! 
H98l IHowes et all 120061: ISchekochihin et all 120090 . For 
Alfvenic fluctuations, the anisotropy implies that, even 
at scales k±pi ~ 1, the turbulent fluctuation frequency 
remains much smaller than the ion cyclotron frequency, 
uj <C fl*. This is an important limit for the applicability 
of gyrokinetic theory, and enables one to determine quan- 
titatively the limit of validity of the heating prescription 
(see iJ2.5l) . When this limit is satisfied, the ion cyclotron 
resonance plays a negligible role in the collisionless d amp- 
ing of the turbulent fluctuations (|Lehe et al.l 120091) . In- 
stead, collisionless damping of the fluctuations occurs via 
the both the ion and electron Landau resonances at scales 
k±p.i > 1, implying that the kinetic Alfven wave cascade 
comprises the "dissipation range" of Alfvenic turbulence. 

It is in the dissipation range that wave-particle inter- 
actions transfer the electromagnetic fluctuation energy 
to the ion and electron particle distribution functions. 
Ultimately, this free energy in the particle distribution 
functions is transferred to sm all scales in velocity space 
through an entropy cascade (ISchekochihin et al.l [20091: 
Tats uno et al 2009; Pl unk et al.ll2010HPhink fc Tatsunol 



201 Hi. enabling arbitrarily weak collisions to thermal- 



ize this energy, increasing the entropy and leading to 
irreversible heating of the plasma. Wave-particle in- 
teractions via the Landau resonance typically peak at 
k±Pi ~ 1 for ions and k±pi > 1 for electrons, so the 
predicted result for ion-to-electron heating Qi / Q e is sen- 
sitive to the model for the nonlinear energy transfer rate 
in the kinetic Alfven wave cascade. 

Although the physical model of the turbulent cascade 
presented here remains controversial within the helio- 
spheric physics community, there exists significant nu- 
merical and observational evidence in support of two of 
its key features: (1) the turbulent frequency remains 
low, uj -C S^i, even for k±pi > 1; and (2) the turbu- 
lence transitions to a cascade of kinetic Alfven waves at 
kj_Pi ~ !• A gyrokinetic numerical simulation of the 
transition from the MHD Alfven wave to the kinetic 
Alfven wave cascade at k±pi > 1 (|Howes et al.l l2008b) 
produces magnetic and electric energy spectra that are 
consistent with Cluster measure ments of turbulence in 
the solar wind (|Bale et al.l [2005V A recent gyrokinetic 
simulation spanning the entire dissip ation range from the 
ion to the electron Larmor radius (Ho wes et al.l 120 111 ) 
yields a magnetic energy spectrum that is quantitatively 
consistent with in situ measureme nts of the dissipation 
range turbulence up to 100 Hz (ISahraoui et al.l 12(3091: 
Kivani et ail 120091: lAlexandrova et al.l 12009b iChen et al.l 



20101: ISahraoui et alJl2010l ). The striking agreement be- 



tween the predictions for a kinetic Alfven wave cascade 
and the o bserved magnet i c and electric power spectra 
found by ISahraoui et al.l ([2009) provide observational 
support for this model. Finally, a fc-filtering analysis 
of multi-spacecraft Cluster measurements demonstrates 



that the wavevectors of the turbulent fluctuations at 
scales k±pi ~ 1 are aligned nearly perpen dicular to the 
local magnetic field (|Sahraoui et al.l [2010T) : for Alfvenic 
turbulent fluctuations, this implies low turbulent fre- 
quencies, uj <C in support of the turbulent model 
employed in this study. 

2.3. Turbulent Heating Prediction 

Assuming a fully ionized plasma of protons and elec- 
trons with i sotropi c Maxw ellian equilibrium velocity dis- 



tributi ons. IHowes! (2010|) used the turbulent cascade 



model (jHowes et al.l 



2008a]) to calculate the total proton 



and electron heating resulting from collisionless damping 
of the electromagnetic fluctuations of the Alfvenic turbu- 
lent cascade. The model empl oys the linear collis ionless 
gyrokinetic dispersion relation ()Howes et al.ll2006l ) to de- 
termine both the linear kinetic damping rate via the Lan- 
dau resonances and the nonlinear energy cascade rate, so 
the resulting heating prescription is only valid in the gy- 
rokinetic limit, uj <C fi p . 

The resulting prescription for the ratio of proton-to- 
electron heating for T p /T e > 1 is given by 



}pl Qe 



Cl 



4 + ^\ m e T, 



rripTp y 



(5) 



where d = 0.92, c 2 = 1.6/(T p /T e ), c 3 = 18 + 
51og(T p /T e ), and a = 2 - 0.2 log(T p /T e ). The require- 
ment that the proton cyclotron resonance plays a negligi- 
ble role in the dynamics and dissipation of the turbulence 
enables the regime of validity of this turbulent heating 
prescription to be quantified. The limit of the regime of 
validity can be expressed as a constraint on the minimum 
width of the inertial range (kop p )^ in , shown in panel (a) 
of Figure[2]as a contour plot in the (/3 p , T p /T e ) parameter 
space. 

Substituting the values of f3 p (R) and T p /T e (R) speci- 
fied in H2.1\ into equation ([5]) enables the calculation of 
the predicted ratio of proton-to-electron turbulent heat- 
ing Q p /Q e {R) as a function of heliocentric radius. In 
Figure [TJ we plot the predicted ratio of proton-to-total 
turbulent heating Q P /{Q P + Q e ) vs. heliocentric radius 
R. 

We compare this theoretical pred i ction to the the em- 
pirical estimate by iCranmer et al.l (|2009f ) (dashed) for 
0.8 AU < R < 5.4 AU in Figure [Q The error esti- 
mates (dotted) in this figure are derived by attempting 
to account for the error arising from both modeling and 
observational uncertainties. Modeling uncertainties are 
derived by taking the curves for outflow speeds of 650, 
700, 750, and 800 km s^ 1 in Figure 3(a) and colatitudes 
of 0, 1 5, and 30 degrees in Figure 4(b) of ICranmer et all 
( 2009) . The observational uncertainties in the calculated 
turbulent heating due to observed variations in the pro- 
ton and electron temperatures and electron heat flux are 
estimated roughly by taking ±10% of the proton-to-total 
turbulent heating ratio. The error estimates (dotted) 
plotted in Figure [1] are determined by taking the outer 
envelope of these modeling and observational uncertain- 
ties. 

We find generally good agreement between the predic- 
tion of the turbulent heating prescription (s olid) and the 
empirical estimate by ICranmer et ahl (|2009f ) (dashed) for 
0.8 AU< R < 5.4 AU. The disagreement at R < 0.8 AU, 
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Fig. 1. — A comparison of the proton-to-total heating rate 
Qp/{Qp + Qe) as a function of heliocentric radius R: the solid 
line is the prediction b ased on the t urbulent heating prescription 
given by equation (Howes 2 0101 ) , and the dashed line is the 
empirical estimate of Cran mer et alj (120091 ') , where dotted lines in- 
dicate an error estimate to account for modeling and observational 
uncertainties. 

as we shall see in H2.5I is attributed to the violation of 
the gyrokinetic approximation, and, therefore, to exceed- 
ing the limits of validity of the turbulent heating pre- 
scription. The downturn in the empirical estimate of 
Qpl (Q P + Qe) (dashed) seen at R > 3 AU may be an arti- 
fact of the bifurcation of the electron tem peratures mea- 
sured by Ulysses, as seen in Figure 1(a) of lCranmer et aLI 
(2009), and may not represent an actual decrease in the 
proto n-to-total turbulent he ating ratio for the high-speed 
wind. ICranmer et al.l (|2009l ) noted that this appeared to 
be a solar cycle effect, but further work will required to 
ascertain the significance of this downturn. 

We can compare our prediction of the turbulent heat- 
ing based on equation (J5J) with predictions based on 
simple theoretical m odels of the turbulent heating by 
ICranmer et al.l (|2009l ). presented in Figure 5 of their pa- 
per. Their approach used a quasilinear framework to 
estimate the proton-to-electron heating rates for three 
particular models of the distribution of turbulent energy 
in wavevector space: an isotropic distribution, a slab 
distribution of only parallel wavevectors, and a "two- 
dimensional" (2D) distribution of nearly perpendicular 
wavevectors. All three of the models showed significant 
disagreement with the empirically determined heating 
ratio. The slab model predicted 100% proton heating, 
the isotropic model overestimated the proton heating for 
R < 2 AU, and the 2D model significantly underesti- 
mated the proton heating at R > 1 AU. In light of these 
results, the agreement between the proton-to-total tur- 
bulent heating ratio predicted by equation ([5]) and the 
empirical estimate in Figure Q] is quite good . This result 
sugges ts that the turbulent cascade model (jHowes et al.1 
2008a) captures the dominant physical mechanisms (de- 
scribed qualitatively in H2.2|) that play a role in the dis- 
sipation of solar wind turbulence at > 0.8 AU. 

The disagreement between the prediction of the 
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6 x 10" 3 


1.2 x 10 5 


5.2 x 10~ 5 
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1.5 x 10" 3 


4.7 x 10 5 


1.3 X 10" 5 


0.9 


6 x 10" 4 


1.2 x 10 6 


5.2 x 10~ 6 



TABLE 1 

Measured values of the frequency of the outer scale 
of the inertial range of solar wind turbulence measured 
BY Helios 2. 



proton-to-total turbulent heating ratio (solid) and the 
empirical estimate (dashed) in Figure [T] can be under- 
stood if we evaluate the limit of the regime of validity 
of the turbulent heating prescription (eq. [5]) using the 
plasma parameters specified for the solar wind model in 
§2.11 The violation of the gyrokinetic approximation, 
or, equivalently, the point at which the proton cyclotron 
resonance begins to play a non-negligible role in the dy- 
namics and dissipation of the turbulence, is cast as a re- 
quirement for the minimum dynamic range spanned by 
the inertial range, given by the driving scale divided by 
the proton Larmor radius, Lq/ p p ~ (fcoPp) We de- 
note this as the minimum "width" of the inertial range, 
{kopp)min> Pitted in panel (a) of Figure [2] as a contour 
plot in the (/3 P , T p /T e ) parameter space. The proton Lar- 
mor radius p p = v tp /fl p = cftTpmp) 1 / 2 /(q p B) is easily 
determined as a function of heliocentric distance R given 
the models for T p (R) and B{R) given in §2.11 Estimating 
isotropic driving wavenumber of the turbulence ko(R), 
however, requires the incorporation of additional empir- 
ical constraints. 

2.4. Evolution of the Width of the Inertial Range 

We interpret the isotropic driving scale Lo = 27r/fco 
to be the outer scale of the inertial range, and we iden- 
tify this scale observationally as the break in the solar 
wind magnetic energy spectrum from the f^ 1 energy 
containing range to the f~ 5 ^ 3 inertial range. To esti- 
mate /co, we employ measurements of magnetic energy 
spectrum from He l ios 2 data published in Figure 23 of 
iBruno fc Carbonel (|2005l ). The frequency /o of the spec- 
tral break marking the outer scale of the inertial range 
measured from this figure is given in Table [TJ Each of 
these spectra are measured in the same corotating fast 
stream with a velocity v sw ~ 7 00 km s _1 (taken from 
Figure 17 of IBruno fc Carbonel ([2005l )K so we may cal- 
culate the corresponding inertial range outer scale length 
Lq = v sw / fo, or wavenumber fco = 2tt/Lq. 

The function 

*=*o(^y (6) 

provides a reasonable fit for the evolution of fco as a 
function of heliocentric radius R with the value K n = 
5 x 10~ 6 rad/km. Although this fit is based solely on 
measurements in the inner heliosphere and may not be 
an accurate representation of the evolution of k (R) for 
R > 0.9 AU, we will see that error in the estimation of 
fco at R > 1 AU does not strongly impact the applicabil- 
ity of the turbulent heating prescription for the plasma 
parameters derived from this solar wind model. 
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Fig. 2. — (a) Logarithmic contour plot of the minimum width of 

the inertial range (kop P )Z^ in over the plane (/3p, T p /T e ), the con- 
dition required for the the proton cyclotron resonance to play a 
negligible role. The red line den otes the evolution plasma param- 
eters of the solar wind model in £|2.1l from R = 0.29 AU (left end) 
to R = 5.4 AU (right end), (b) The minimum width of the inertial 
range (kopp)~ • along the red line in panel (a) (dashed) vs. the 
estimated width of the inertial range (kop p )~ 1 for the solar wind 
model considered here (solid). For the heating prescription to be 
valid, the dashed line must fall below the solid line. 

2.5. Evaluation of the Limits of Validity of the 
Turbulent Heating Prescription 

Using the values of p p (R) and ko(R) derived here, 
we calculate the width of the inertial range k$p p as a 
function of heliocentric radius and compare it to the 
constraint on {kap p )~^ in , as shown in Figure [2] Panel 
(a) presents a logarithmic contour plot of the constraint 

(fcoPp)min over t ne (Pp,T p /T e ) parameter space. Also 
shown is the path though this parameter space (red) 
traversed by the solar wind model from 0.29 AU to 
5.4 AU. In panel (b), the value of the width of the in- 
ertial range (fc p p ) _1 for our solar wind model (solid) is 
plotted against the constraint on the minimum width of 
the inertial range (kop P )^ in for the validity of our heat- 
ing model (dashed). (For the heating prescription to be 
valid, the dashed line must fall below the solid line in 
this plot.) As previously mentioned, for the plasma pa- 
rameters P p and T p /T e specified by the solar wind model 

in £12.11 the constraint on (fcoPp)min does not strongly 
restrict the applicability of the heating prescription for 
R > 1 AU. It is clear, however, that this constraint is 
violated for heliocentric distances R < 0.7 AU, signaling 
that the gyrokinetic representation of the turbulent dissi- 
pation mechanisms is no longer valid because the proton 
cyclotron resonance has begun to play a non-negligible 
role. Significantly, this limit to the validity of our predic- 
tion for the proton-to-total turbulent heating ratio coin- 
cides with the point where our prediction (solid) begins 
to fall significantly below the empirical estimate (dashed) 
in Figure [T] If the contribution to proton heating from 
the proton cyclotron resonance becomes non-negligible 
for heliocentric distances R < 0.8 AU, a physical pro- 



cess not represented in the gyrokinetic turbulent heat- 
ing prescription, then the result would be an empirically 
estimated proton-to-total heating rate that exceeds the 
gyrokinetic predictions at these radii, in agreement with 
the behavior shown in Figure [T] 

3. DISCUSSION 

The results presented in Ej2] suggest the following con- 
sistent picture of the physical mechanisms in the high- 
speed solar wind that are responsible for the dissipation 
of the turbulence and that lead to heating of the plasma 
protons and electrons. 

In the inner heliosphere at R < 0.8 AU, the typi- 
cally high T p /T e and low f3 p conditions (coupled with 
a slightly smaller width of the inertial range (fcoPp) 1 ; 
see panel (b) of Figure ^ in high-speed streams lead to 
a turbulent cascade in which the small scale turbulent 
fluctuations can reach the proton cyclotron frequency 
before the turbulence is collisionlessly damped via the 
Landau resonances. Therefore, one may expect to ob- 
serve greater heating of the protons than that predicted 
by the turbulent heating prescription given by equa- 
tion (Et, a result ba sed on a gyrokinetic cascade model 
(jHowes et al.l I2008al ) . By the time the turbulence has 
reached R > 0.8 AU, the decrease of T p /T e and increase 
of f3 p (coupled with a slight increase in the width of 
the inertial range (fcoPp) -1 ; see panel (b) of Figure [2]), 
lead to plasma conditions in which the turbulent cas- 
cade no longer is affected by the proton cyclotron res- 
onance before it is terminated by collisionless damping 
via the Landau resonances. Thus, for the range 0.8 AU 
< R < 5.4 AU, the turbulent heating prescription in the 
gyrokinetic limit adequately represents all of the physi- 
cal mechanisms needed to reproduce the empirically es- 
timated proton-to-total turbulent heating ratio, as seen 
in Figure [T] 

This theoretical prediction of non-negligible proton cy- 
clotron damping within the inner heliosphere at R < 
0.8 AU is consistent with observational evidence. Proton 
cyclotron damping is expected to lead to heating of the 
protons in the direction perpend icular to the local mean 
magnetic field (jLehe et "all [20 09). In the absence of such 
perpendicular proton heating, double adiabatic evolution 
would lead to a const ant value of T± v / B as a function 
of heliocentric radius (|Chew et al.l [l956). Helios obser- 
vations demonstrate that the constancy o f Tj n/B is in- 
deed violated within the inner heliosphere (jMarsch et al.l 
[198! . 

It is important to note that the observed non-adiabatic 
T± p /B does not directly identify the physical mecha- 
nism responsible. In addition to heating via the pro- 
ton cyclotron resonance, several other physical mech- 
anisms could lead to the observed behavior, including 
stochastic proton heating and kinetic proton tempera- 
ture anisotropy instabilities. 

Stochastic heating had been proposed as a mecha- 
nism for perpendicular proton heating for some time 
(I Johnson fc Chend [20011 iChen et all 12001 1 White et all 
2002; Voiten ko fc Gpossensl l2004t iBourouaine et aLl 
120081 ). and recently iChandran et al.l ([ToiOD have put 
forth strong theoretical and numerical evidence for 
stochastic perpendicular heating of protons at low 
plasma by kine tic Alfven waves of sufficient amplitude. 
iChandranl (|2010P ) has used these results for stochastic 
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heating to explain proton and minor ion perpendicular 
temperature observations in coronal holes. For the low 
plasma j3 conditions found in the inner heliosphere, it is 
possible that this mechanism of stochastic heating could 
explai n the perpendicular proton temperature measure- 
ments (|Marsch et al . 1983) and the em pirically estimated 
proton heating (jCranmer et al.ll2009D at R < 0.8 AU. 

In addition, kinetic proton temperature anisotropy 
instabilities in the spherically expanding solar wind 
flow have been shown to play a role in regulating 
the proton temperat ure anisotro py (|Kasper et al.l [20021 : 
Bellinger et"all l2006t iBale et al.l 120091 ) . Although these 
instabilities cannot lead to a net heating of the pro- 
ton species, they can mediate a transfer of energy from 
the parallel to the perpendicular temperature, and vice 
versa, and so may be a cause of the observed devia- 
tion from double adiabatic e volution of the perpe ndicular 
temperature of the protons ([Marsch et aljfl983l ). 

There are several possible limitations of the applica- 
tion of the turbulent heating prescription given by equa- 
tion ([5]) to the problem of the heating due to the dissi- 
pation of solar wind turbulence. First, the model does 
not account for energy in compressible wave modes, such 
as the collisionless manifestation of the fast and slow 
MHD wave modes. If significant energy exists in these 
compressible modes, any heating due to the dissipa- 
tion of these modes must be handled separately. Sec- 
ond, the cascade model is constructed specifically for 
the case of balanced Alfven wave energy fluxes up and 
down the local magnetic field, corresponding to zero 
cross helicity. Observations of the cross helicity in the 
high-latitude wind typically show non-zero normalized 
cross helicities varying over the range 0.2 < a c < 0.6 
(jBavassano et al.ll2000al |bT). Note, however, that a cross 
helicity of o~ c = 0.6 corresponds to the amplitude of anti- 
sunward waves only a factor of 2 larger than the sun- 
ward waves. Since the typical energy cascade rates in 
strong M HD turbulence vary linear l y with the wave am - 
plitudes (|Goldreich fc SridhariU9§5t faowes et al.ll2008al) . 
this level of imbalance in the turbulence is unlikely to 
yield significant qualitative differences compared to the 
balanced case. Third, the simplifying assumption of a 
fully ionized proton and electron plasma with isotropic 
Maxwellian equilibrium velocity distributions neglects 
the physical variations that may arise from the more 
complicated equilibrium conditions often observed in 
the solar wind. Such conditions include temperature 
anisotropy with respect to the local magnetic field di- 
rection (often treated using a bi-Maxwellian equilibrium 
distribution), significant deviations from a Maxwellian 
distribution at high energy, and the presence of minor 
ions, particularly helium. Nonetheless, we believe that 
the results presented here represent a significant step for- 
ward in our understanding of the mechanisms responsi- 
ble for proton and electron heating in the turbulent solar 



wind. Future work will explore the implications of these 
additional effects if any of them appears to impact sig- 
nificantly the findings presented here. 

4. CONCLUSIONS 

In the effort to identify the physical mechanisms that 
govern the dissipation of solar wind turbulence and 
lead to heating of the solar wind protons and electrons, 
iCranmer et al.l (|2009f ) made a great stride forward by 
determining an empirical estimate of the proton-to-total 
plasma heating in the high-speed solar wind using He- 
lios and Ulysses data. Based on a turbulent energy cas- 
cade model for low-frequency, anisot ropic Alfvenic turbu - 
lence in a wea kly collisional plasma (|Howes et al.l l2~008a) , 
iHowesi (|2010l ) constructed an analytical prescription for 
the total proton-to-electron heating resulting from colli- 
sionless damping of the electromagnetic fluctuations of 
the Alfvenic turbulence. 

Applying this turbulent heating prescription to predict 
the proton-to-total plasma heating in the high-speed so- 
lar wind, we obtain the following results, as shown in 
Figure [TJ (1) the prediction agrees well with the empir- 
ical estimate for 0.8 AU < R < 5.4 AU; (2) the pre- 
dicted proton heating falls below the empirical estimate 
for R < 0.8 AU. Investigating the cause of the disagree- 
ment for R < 0.8 AU, we see, in panel (b) of Figure [21 
that the turbulent heating prescription in the gyrokinetic 
limit ceases to be valid for R < 0.7 AU. This failure of the 
prescription's validity has physical meaning. In this re- 
gion, the cascade model predicts that the turbulent fluc- 
tuations will reach the proton cyclotron frequency before 
they are damped via the Landau resonances. Therefore, 
we expect that proton cyclotron damping will cause ad- 
ditional proton heating, leading to an underestimate of 
the proton-to-total heating ratio, as seen in Figure [T] 

These results suggest the following physical picture of 
the turbulent cascade and plasma heating in the high- 
speed solar wind. In the inner heliosphere at R < 0.8 AU, 
the turbulent cascade can reach the proton cyclotron 
frequency, leading to a level of proton heating higher 
than that predicted by the turbulent heating prescrip- 
tion in the gyrokinetic limit. But for heliocentric radii 
R > 0.8 AU, collisionless damping via the Landau reso- 
nances terminates the turbulent cascade before the pro- 
ton cyclotron resonance is reached, so the turbulent heat- 
ing prescription in the gyrokinetic limit contains all of the 
necessary physical mechanisms needed to reproduce the 
empirically estimated proton-to-total heating ratio. 

G. G. H. thanks Steve Cranmer for insightful dis- 
cussions and for providing his empirical heating data 
with error estimates. The work has been supported 
by NSF CAREER Award AGS-1054061 and NASA 
NNX10AC91G. 



REFERENCES 



Alexandrova, O., Saur, J., Lacombe, O, Mangeney, A., Mitchell, 
J., Schwartz, S. J., & Robert, P. 2009, Phys. Rev. Lett., 103, 
165003 

Bale, S. D., Kasper, J. O, Howes, G. C, Quataert, E., Salem, C, 
& Sundkvist, D. 2009, Phys. Rev. Lett., 103, 211101 

Bale, S. D., Kellogg, P. J., Mozer, F. S., Horbury, T. S., & Reme, 
H. 2005, Phys. Rev. Lett., 94, 215002 



Bavassano, B., Pietropaolo, E., & Bruno, R. 2000a, 

J. Geophys. Res., 105, 12697 
— . 2000b, J. Geophys. Res., 105, 15959 
Boldyrev, S. 2005, Astrophys. J. Lett., 626, L37 
Bourouaine, S., Marsch, E., & Vocks, C. 2008, 

Astrophys. J. Lett., 684, L119 



Solar Wind Turbulent Heating 



7 



Breech, B., Matthaeus, W. H., Cranmer, S. R., Kasper, J. C, & 
Oughton, S. 2009, Journal of Geophysical Research (Space 
Physics), 114, 9103 

Breech, B., Matthaeus, W. H., Minnie, J., Bieber, J. W., 

Oughton, S., Smith, C. W., & Isenberg, P. A. 2008, Journal of 
Geophysical Research (Space Physics), 113, 8105 

Bruno, R., & Carbone, V. 2005, Living Reviews in Solar Physics, 
2, 4 

Chandran, B. D. G. 2010, Astrophys. J., 720, 548 
Chandran, B. D. G., Li, B., Rogers, B. N., Quataert, E., & 

Germaschewski, K. 2010, Astrophys. J., 720, 503 
Chen, C. H. K., Horbury, T. S., Schekochihin, A. A., Wicks, 

R. T., Alexandrova, O., & Mitchell, J. 2010, Physical Review 

Letters, 104, 255002 
Chen, L., Lin, Z., & White, R. 2001, Phys. Plasmas, 8, 4713 
Chew, G. L., Goldberger, M. L., & Low, F. E. 1956, 

"Proc. R. Soc. London A", 236, 112 
Cranmer, S. R., Matthaeus, W. H., Breech, B. A., & Kasper, 

J. C. 2009, Astrophys. J., 702, 1604 
Frieman, E. A., & Chen, L. 1982, Phys. Fluids, 25, 502 
Goldreich, P., & Sridhar, S. 1995, Astrophys. J., 438, 763 
Gruzinov, A. V. 1998, Astrophys. J., 501, 787 

Hellinger, P., Travm'cek, P., Kasper, J. C, & Lazarus, A. J. 2006, 

Geophys. Res. Lett., 33, 9101 
Howes, G. G. 2008, Phys. Plasmas, 15, 055904 
— . 2010, Mon. Not. Roy. Astron. Soc, 409, L104 
Howes, G. G., Cowley, S. C, Dorland, W., Hammett, G. W., 

Quataert, E., & Schekochihin, A. A. 2006, Astrophys. J., 651, 

590 

— . 2008a, J. Geophys. Res., 113, A05103 

Howes, G. G., Dorland, W., Cowley, S. C, Hammett, G. W., 

Quataert, E., Schekochihin, A. A., & Tatsuno, T. 2008b, 

Phys. Rev. Lett., 100, 065004 
Howes, G. G., TenBarge, J. M., Dorland, W., Quataert, E., 

Schekochihin, A. A., Numata, R., & Tatsuno, T. 2011, Phys. 

Rev. Lett., submitted 
Johnson, J. R., & Cheng, C. Z. 2001, Geophys. Res. Lett., 28, 

4421 

Kasper, J. C, Lazarus, A. J., & Gary, S. P. 2002, 
Geophys. Res. Lett., 29, 20 



Kiyani, K. H., Chapman, S. C, Khotyaintsev, Y. V., Dunlop, 

M. W., & Sahraoui, F. 2009, Phys. Rev. Lett., 103, 075006 
Kolmogorov, A. N. 1941, Dokl. Akad. Nauk SSSR, 30, 9, english 

Translation: Proc. Roy. Soc. London A, 434, 9 (1991) 
Leamon, R. J., Smith, C. W., Ness, N. F., & Wong, H. K. 1999, 

J. Geophys. Res., 104, 22331 
Lchc, R., Parrish, I. J., & Quataert, E. 2009, Astrophys. J., 707, 

404 

Marsch, E., Muchlhaeuser, K. H., Rosenbauer, H., & Schwenn, R. 

1983, J. Geophys. Res., 88, 2982 
Matthaeus, W. H., Zank, G. P., Smith, C. W., & Oughton, S. 

1999, Physical Review Letters, 82, 3444 
Plunk, G. G., Cowley, S. C, Schekochihin, A. A., & Tatsuno, T. 

2010, Journal of Fluid Mechanics, 664, 407 
Plunk, G. G., & Tatsuno, T. 2011, Phys. Rev. Lett., 106, 165003 
Quataert, E. 1998, Astrophys. J., 500, 978 
Quataert, E., & Gruzinov, A. 1999, Astrophys. J., 520, 248 
Rutherford, P. H., & Frieman, E. A. 1968, Phys. Fluids, 11, 569 
Sahraoui, F., Goldstein, M. L., Belmont, G., Canu, P., & Rezeau, 

L. 2010, Phys. Rev. Lett., 105, 131101 
Sahraoui, F., Goldstein, M. L., Robert, P., & Khotyaintsev, Y. V. 

2009, Phys. Rev. Lett., 102, 231102 
Schekochihin, A. A., Cowley, S. C, Dorland, W., Hammett, 

G. W., Howes, G. G., Quataert, E., & Tatsuno, T. 2009, 

Astrophys. J. Supp., 182, 310 
Smith, C. W., Matthaeus, W. H., Zank, G. P., Ness, N. F., 

Oughton, S., & Richardson, J. D. 2001, J. Geophys. Res., 106, 

8253 

Stawarz, J. E., Smith, C. W., Vasquez, B. J., Forman, M. A., & 

MacBride, B. T. 2009, Astrophys. J., 697, 1119 
Tatsuno, T., Schekochihin, A. A., Dorland, W., Plunk, G., 

Barnes, M. A., Cowley, S. C, & Howes, G. G. 2009, 

Phys. Rev. Lett., 103, 015003 
Voitcnko, Y., & Goossens, M. 2004, Nonlinear Processes in 

Geophysics, 11, 535 
White, R., Chen, L., & Lin, Z. 2002, Phys. Plasmas, 9, 1890 
Zank, G. P., Matthaeus, W. H., & Smith, C. W. 1996, 

J. Geophys. Res., 101, 17093 



